Table 6. Spectral fitting results for NGC3690. 

Parameter Power-law Single Temperature Double Temperature Raymond Smith + Triple temperature ce\ 

Raymond-Smith Raymond Smith Power-law Raymond Smith 

kT(KeV) 4.25±°'^ 0.83±°-°| 0.83±°;°4 57.05 

10.3t£S 0.871+°;^ 

r 9 9 o + 0.12 , r- fi +0.11 

1 Z.Za_ Q 12 J - OD — 0.11 

at 0.22 

7V H (10 2 ° C m- 2 ) 5.50+0-; 2.11+°;|| l-60±g;g 2.42+°;« 2.031+°; 2 ^ 3.2] 

X 2 / d.o.f. 490.0/242 624.5/242 287.8/235 291.5/235 290.3/237 309. 

Flux* (0.1-2.0keV) 9.72 8.85 8.5 9.0 8.95 8 

Flux* (2.0-10.0keV) 7.68 9.59 10.8 10.8 10.3 ! 

Total Luminosity* (0.1-10.0keV) 8.3 8.8 9.2 9.46 9.2 

QQ i f The slope of power-law distribution of temperatures in the cevmekl model. 

■ | The fluxes are in units of 10 _13 ergs — 1 cm -2 . 

' * The luminosity is in units of 10 41 ergs — 1 . 
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ABSTRACT 

We present ROSAT (HRI and PSPC) and ASCA observations of the two luminous 



10 4 



erg s" 1 ) star-forming galaxies NGC3310 and NGC3690. The HRI 



shows clearly that the sources are extended with the X-ray emission in NGC3690 
coming from at least three regions. The combined 0.1-10 keV spectrum of NGC3310 
can be described by two components, a Raymond-Smith plasma with temperature 



kT = 0.81_qi2 keV and a hard power-law, T = 1.44_o' n , (or alternatively a harder 
Raymond-Smith plasma with kT ~ 15 keV), while there is no substantial excess 
absorption above the Galactic. The soft component emission is probably due to a 
super- wind while the nature of the hard emission is more uncertain with likely origins, 
X-ray binaries, inverse Compton scattering of IR photons, an AGN or a very hot gas 

+0.02 
04 



component (~ 10 8 K). The spectrum of NGC3690 is similar, with kT = 0.83Z V 



keV and T = 1.56 



+0.11 
-0.11' 



We also employ more complicated models such as a multi- 



temperature thermal plasma, a non-equilibrium ionization code or the addition of 
a third softer component which improve the fit but not at a statistically significant 
level (< 2a). These results are similar to recent results on the archetypal star- forming 
galaxies M82 and NGC253. 

Key words: galaxies: starburst-galaxies: merger-X-rays: galaxies-galaxies: individual: 
NGC3690-galaxies: individual: NGC3310 



1 INTRODUCTION 

Significant X-ray emission from star-forming galaxies has 
been known since the EINSTEIN epoch (eg Stewart et al. 
1982). Their X-ray spectrum in the 0.3-3.5 keV band can be 
fitted by a thermal plasma of temperature about 2 keV with 
an absorption often higher than the Galactic value (eg. Fab- 
biano 1988, Kim et al. 1992). However, simple spectral fits 
suggested that the emission mechanism may be more com- 
plicated than a single temperature thermal plasma. Indeed, 
for a few nearby star-forming galaxies (eg M82 and NGC253) 
several discrete sources, mostly X-ray binaries, are also re- 
solved (Long & Van Speybroeck 1983). The EINSTEIN HRI 
provided the first evidence for extended X-ray features (eg 
Watson et al. ,1984, Fabbiano and Trinchieri 1984). 

The effective area of ROSAT PSPC provided the op- 
portunity to dramatically increase the number of starform- 
ing galaxies studied in X-rays. Hot gas halos and outflows, 
have also been observed, around many starforming galaxies 
(eg. Read, Ponman & Wolstencroft 1995, Delia Ceca, Grif- 
fiths & Heckman 1996, Read, Ponman & Strickland 1997), 
as was anticipated based on previous optical (Heckman et 
al. 1990) and theoretical work (e.g. MacLow and Mc Cray 



1988). This wind can be produced by the supply of mechan- 
ical energy via stellar winds from evolved massive stars and 
numerous supernovae. The hot gas is accelerated outwards 
forming an expanding superbubble, which on reaching the 
extent of the minor axis of the galaxy, blows out following 
the onset of the Rayleigh- Taylor instability. These outflows 
can extend for several kpc along the minor axis of the galaxy, 
and have major influence on its subsequent evolution. This 
is especially so in the case of dwarf galaxies, where because 
of their weak gravitational potential, the wind can remove 
the interstellar medium causing star formation to cease, (eg 
Dekel & Silk 1986, Heckman et al. 1995). Combined ROSAT 
and ASCA observations of star-forming galaxies, (eg Ser- 
lemitsos, Ptak & Yaqoob 1996) suggest that the soft X-ray 
band can be fitted with a thermal plasma and a power-law 
model. The plasma temperature is ~ 0.6 — 0.8keV, with an 
absorbing column Nh > 10 21 cm -2 , well above the Galactic 
value, and the power-law photon index is ~ 1.7. The hard X- 
ray band can be fitted equally well with a high temperature 
thermal plasma model, kT > 6keV, (cf. Ohashi et al. 1990). 
The soft emission may emanate from the super-wind while 
the origin of the harder power-law X-ray emission may be 
due to either X-ray binaries or Inverse Compton scattering 
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of infrared-photons by the relativistic electrons generated 
by the supernovae (eg Rieke et al. 1980). The presence of 
very hot gas (~ 10 8 K), or a weak active nucleus cannot be 
ruled out. However, medium spectral resolution, wide energy 
band (0.5-10 keV) observations of NGC253 and M82 with 
ASCA, show that the above two component models do not 
give an acceptable description of the data (Ptak et al 1997, 
Moran & Lehnert 1997). Even the introduction of a third 
thermal component (kT ~ 0.3keV) still provides a poor fit to 
the high-resolution ASCA SIS spectrum (Moran & Lehnert 

1997) . This suggests that the emission mechanisms are more 
complex than these simple models can explain. Since the X- 
ray spectra of Galactic supernova remnants (eg Hughes & 
Singh 1994) cannot be well-fit by a simple Raymond-Smith 
model (their emission is dominated by non-equilibrium pro- 
cesses) it is likely that we are dealing with a similar situation 
in starforming galaxies. 

One of the interesting results of the ROSAT All-Sky 
Survey was the discovery of star-forming galaxies with X- 
ray luminosities ~ 10 42 ergs -1 , approaching those of Seyfert 
galaxies (eg Moran et al 1995). In order to gain further in- 
sight into the nature of this powerful X-ray emission and 
to compare starburst properties with those of the archety- 
pal (less luminous) star-forming galaxies, we have cross- 
correlated the ROSAT All Sky Survey Bright Source Cata- 
logue, RASSBSC, (Voges et al. 1996) with the spectral atlas 
of nearby galaxies of Ho, Filippenko and Sargent (1995). The 
advantage of the latter sample is that it has high quality op- 
tical spectra and thus provides us with a reliable classifica- 
tion for each galaxy. The sample of Ho et al. contains galax- 
ies with declination S > 0° and magnitude B < 12.5 while 
the RASS contains almost 20,000 sources over the whole sky. 
The cross-correlation yielded 43 X-ray counterparts within 
1 arcmin of the optical galaxy, whereas less than one coinci- 
dence is expected by chance. Within this sample of galaxies 
the vast majority are AGNs (Seyferts 1&2, and LINERS), 
but there are also some star-forming galaxies and high lumi- 
nosity early-type galaxies. Further details of the properties 
of the resulting X-ray sample are given elsewhere (Roberts, 

1998) . The star- forming galaxies are listed in table 1. In this 
paper we present X-ray imaging and spectral analysis for 
the two most luminous star-forming galaxies in the sample: 
NGC3310 and NGC3690, excluding NGC5905 which under- 
went some form of transient phenomenon as yet not under- 
stood. Both ROSAT and ASCA observations are available 
in the HEASARC archive for these two galaxies. 

We use a value of H = 50km s -1 Mpc -1 for the Hubble 
constant, throughout this paper. All errors quoted refer to 
the 90 per cent confidence level. 

1.1 NGC 3310 

NGC3310 is a well studied nearby galaxy of Sbc(r)pec type 
(Mulder et al. 1996). Its recession velocity is 980km s -1 (dc 
Vaucouleurs et al. 1991) which implies a distance of 19.6 
Mpc. Its optical image shows some interesting features. The 
most prominent is an arc and bow structure which extends 
for 100" to the North West of the centre of the galaxy. 
Bertolla and Sharp (1984) propose that the arc could be part 
of a spiral arm and the bow is the remnant of an old jet. 
Balick and Heckman (1985) suggest that the whole struc- 
ture is the remnant of a collision between the galaxy and a 



Table 1. The X-ray luminous starburst galaxies 
in the Ho et al. sample 



Galaxy X — ray luminosity* 

M82 40.88 

NGC3310 40.89 

NGC3690 41.62 

NGC4449 39.24 

NGC5204 39.54 

NGC5905 42.18* 



f This is the log of the X-ray luminosity in the 
ROSAT band, from the RASS data. 
| This galaxy has show significant long time vari- 
ability (Bade et al. 1996); This luminosity refers 
at the outburst epoch. 



dwarf companion. This is in agreement with the suggestion 
of Mulder et al. (1985) who characterize the system as a 
young merger based on the anomalies found in its rotation 
curve. It should be noted that in the radio and near-infrared 
images we do not see two resolved nuclei (Telesco and Gat- 
ley 1984, Balick and Heckman 1981). This could mean that 
the merger phase is now complete as suggested by Balick 
and Heckman. 

Another interesting morphological feature is a giant HII 
region situated ~ 12" SW of the nucleus (Balick and Heck- 
man 1981). Its size is comparable to the largest extragalactic 
HII regions known, and its spectral characteristics are dom- 
inated by signatures of Wolf Rayet stars (Pastoriza et al. 
1993). 

The starburst nature of the activity has been confirmed 
by numerous observations in many spectral bands (e.g. Tele- 
sco and Gatley 1984, Smith et al. 1996, Balick and Heckman 
1981). Evolutionary synthesis models of the starburst give 
an estimated age range of 10 7 to 10 s years, (Van der Kruit 
and de Bruyn 1976, Balick and Heckman 1981, Telesco and 
Gatley 1984, Pastoriza et al. 1993). 

NGC3310 has been previously observed with the Imag- 
ing Proportional Counter (IPC) on board the EINSTEIN 
observatory (Fabianno et al. 1992) yielding an X-ray flux 
f x = 1.1 x 10 -12 erg s -1 cm -2 in the 0.3-3.5 keV band 
for a 5keV bremsstrahlung model and Galactic absorp- 
tion. This corresponds to an X-ray luminosity of 5.0 x 
10 40 ergs -1 cm" 2 . 

1.2 NGC 3690 

NGC3690 is also a merging system, but in contrast to 
NGC3310 it is in an earlier merger phase. The system is 
composed of NGC3690, which is the western part of the 
merger and IC694, the eastern component. The whole sys- 
tem is also known as Arp 299 or Mrk 171. Casoli et al. (1989) 
propose that there is a third component in the system thus 
forming an interacting triplet. Its recession velocity is 3159 
km s -1 (Sanders and Mirabel 1985), implying a distance of 
63.2 Mpc. At that distance the projected separation of 22.5" 
of the two components corresponds to 6.9 kpc. The infrared 
surface brightness of NGC3690 is at least twice that of IC694 
(Friedman et al. 1987). The total mass of gas calculated from 
radio observations is about 2 x 10 11 M Q (Casoli et al. . 1989) 
The high IR luminosity (L FIR = 1.2 x 1O 12 L ) (Soifer et 
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al. 1987), results from re-radiating dust heated by the star- 
burst activity. This is supported by multi-wavelength stud- 
ies (Gehrz et al. 1983, Friedman et al. 1987, Nakagawa et al. 
1989). Detailed analysis of emission line ratios using the di- 
agnostic diagrams of Veilleux and Osterbrock (1987) clearly 
indicate a stellar origin of the ionizing continuum (Friedman 
et al. 1987). The infrared emission is extended over a region 
of several kpc, but the main sources are located close to the 
two nuclei (Friedman et al. 1987). Modeling of the starburst 
has been carried out by Gehrz et al. (1983) and Nakagawa et 
al. (1989). The latter find a starburst age of about lOMyrs. 
An interesting point is that the two nuclei have different 
properties, implying either a different age and/or a different 
Initial Mass Function (IMF) for each starburst event. 

NGC3690 has been previously observed with the High 
Resolution Imager (HRI) on board the EINSTEIN obser- 
vatory. 17 counts were detected implying an X-ray flux 
f x = 4.8 x 10 _13 crgs _1 cm -2 in the 0.2-4.0keV band, for a 5 
keV bremsstrahlung model and Galactic absorption (Fabi- 
anno et al 1992). This corresponds to an X-ray luminosity 
of 2.3 x 10 41 ergs^cm- 2 . 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 The ROSAT PSPC Observations 

NGC3310 and NGC3690 were each observed on two oc- 
casions with the Position Sensitive Proportional Counter 
(PSPC) (Pfefferman et al. 1987) on board ROSAT (Triimper 
et al. 1984). The details of each observation are given in ta- 
ble 2. 

For the data reduction we have followed the standard 
procedure, using the ASTERIX package. We excluded those 
data with Master Veto rates higher than 170cps. We ex- 
tracted a PSPC spectral image (channels 11 to 201) with a 
pixel size of 15". To obtain an X-ray spectrum we extracted 
data from a circular region of 1.5' radius around the X-ray 
centroid. The background was estimated from an annular 
region between radii of 15' and 8.8' from the centroid, after 
exclusion of all the discrete sources detected with the PSS 
algorithm (Allen 1992) down to the 5a level. 

2.2 The ROSAT HRI Observations 

Both galaxies have been observed with the High Resolution 
Imager (HRI) (David et al. 1997) on board ROSAT . HRI 
covers a field of 38' diameter and it has a spatial resolution 
of about 5" (the FWHM of the XRT+HRI Point Spread 
Function). We note that the HRI has very limited spectral 
resolution so we cannot use it for spectral analysis (David 
et al. 1997). The screening of the data has been carried out 
using the ASTERIX package. We have rejected all data with 
aspect errors greater than 2. 

2.3 The ASCA Observations 

Both galaxies have been observed with the ASCA satellite 
(Tanaka et al. 1994). On board ASCA there are four instru- 
ments: two Gas Imaging Spectrometers (GIS2 and GIS3) 
(Ohashi et al. 1996) and two Solid-State Imaging Spectrom- 
eters (SISO and SIS1) (Gendreau 1995). For screening of the 
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Figure 1. The HRI radial profile of NGC3310 (circles) overlaid 
on the radial profile of AR Lac (stars). 

data we have followed the standard procedure and used the 
ASCACLEAN program in the FTOOLS package with the 
parameters described in the ABC ASCA Guide (Yaqoob et 
al. 1997). After the standard processing we inspected the 
light curves and removed time intervals with unusually high 
numbers of counts, which may result from background par- 
ticle contamination. The extraction of the source and the 
background spectra was been carried out using the XSE- 
LECT program within the FTOOLS package. We extract 
the source spectrum using a circular region of 2.7 and 5.5 
arcmin radius for the SIS and GIS respectively. For the back- 
ground regions we have selected ten source free zones, each 
of area about 6 and 10 arcmin 2 for the SIS and GIS respec- 
tively. Following extraction of the source and background 
spectra, we ran the SISrmg program to create the SIS re- 
sponse matrices for each SIS source spectral file. For the GIS 
spectral files we used the gis2v4_0.rmf and gis3v4_0.rmf for 
GIS2 and GIS3 respectively. We created the Ancillary Re- 
sponse Matrices for each file by running the ASCAarf pro- 
grams within the FTOOLS package. 



3 SPATIAL EXTENT ANALYSIS 
3.1 NGC3310 

We used the archival HRI observations to obtain informa- 
tion on the morphology of the X-ray emitting regions and 
compared these with images of these galaxies at other wave- 
lengths. First we searched for extended emission by compar- 
ing the background subtracted radial profile of the galaxies 
with the radial profile of a point source. For this purpose 
we used as a reference the point source corresponding to the 
star AR-Lac. We retrived archival data for an on-axis obser- 
vation of AR-Lac, and applied the same extraction process 
as used for the galaxy. The comparison of the radial profiles 
of NGC3310, and the radial profile of AR-Lac is shown in 
figure 1. This figure clearly shows that the X-ray source has 
a radius of ~ 0.5 arcmin (at larger radii the emission from 
the galaxy begins to blend with the profile of AR-Lac, within 
the 2a error bars). Although this is a somewhat arbitrary 
point at which to define the extent of the emission, it gives 
a lower limit to the source size. At the distance of NGC3310 
this corresponds to ~ 3 kpc. We note that fitting a Gaus- 
sian function to the radial profile gives a similar source size 
(FWHM=1.6 arcmin). 
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Table 2. Summary of the observations 
Galaxy 

NGC3310 



NGC3690 



oat-ciiit-c 


Instrument 


Date Observation started 


Exposure Time| 








Ksec 




pepp 
ir or O 


±yy i- 1 1- id 






HRI 


1 QQ^l 04. 17 


41.842 


ROSAT 


HRI 


1994-11-23 


4095 


ASCA 




1994-04-17 


GIS: 10.4, SIS:10.4 


ASCA 




1994-11-13 


GIS:11.6, SIS:10.2 


ROSAT 


PSPC 


1993-04-22 


3.534 


ROSAT 


PSPC 


1991-11-18 


6.391 


ROSAT 


HRI 


1993-04-18 


6.751 


ASCA 




1994-04-06 


GIS:5.7, SIS:5.3 


ASCA 




1994-12-01 


GIS:38.2, SIS:35.5 



In order to compare the distribution of the X-ray emit- 
ting gas with the optical image we have resampled the X-ray 
image with a pixel size of 1.5 arc-seconds and then smoothed 
this image using a two-dimensional gaussian of FWHM=3.5 
arc-seconds (2.3 pixels) following Delia Ceca et al. (1996). 
The final pixel size is 5.4", equal to the XRT+HRI PSF. 
Then we obtained contours from this image corresponding to 
2.2, 2.0 ,1.7, 1.33, 0.88, 0.44, 0.35, 0.26, 0.17 counts arcsec" 2 
and overlaid these onto an optical image obtained from the 
Digitized Sky Survey. This image is shown in figure 2. The 
main problem we faced in attempting to overlay the two im- 
ages is that of frame registration. There are no X-ray sources 
within the HRI image which have an obvious optical coun- 
terpart. Thus in order to overlay the X-ray contour plot onto 
the optical image, we are forced to assume that the centroid 
of the X-ray source corresponds to that of the optical nu- 
cleus. In the X-ray image of NGC3310 we detect two other 
sources in addition to the nucleus, at a significance above 
10a. The details concerning these sources are presented in 
table 3. These point sources may correspond to luminous 
X-ray binaries or young supernova remnants. Unfortunately 
there are no identified optical counterparts for these X-ray 
sources on the POSS plates. 

3.2 NGC3690 

The HRI analysis of NGC3690 followed the same procedure 
as for NGC3310. The image pixels were binned to a size of 
1.5" and then smoothed using a two-dimensional gaussian 
of FWHM 3.5". Finally, we overlaid the X-ray contours on 
a POSS plate image obtained from the Digitized Sky Sur- 
vey, figure 3; the contours correspond to 0.57, 0.44, 0.40, 
0.35, 0.31, 0.27, 0.22, 0.18, 0.13 and 0.11 counts arcsec" 1 re- 
spectively. The most striking feature is the existence of three 
separate sources in the X-ray image. The two most luminous 
sources correspond to the two merging nuclei. Their X-ray 
fluxes are given in table 4. The third X-ray source corre- 
sponds to a hot spot seen in the infrared and radio images 
to the north of NGC3690 (the western component). Apart 
from these three sources there are marginal detections of 
other sources, but at a low level of significance (below 3a), 
and will not be discussed further. We note that the relative 
strengths of the three main X-ray sources follow quite well 
the relative strengths of the near infrared sources, suggesting 
a possible common origin for the IR and the X-ray emission. 

We have searched for extended emission from the two 
resolved nuclei. From figure 3 it is clear that the emission 



from the two resolved nuclei is extended and not symmet- 
rical. Fitting a Gaussian function to the radial profiles of 
NGC3690 and IC694 yields a FWHM of 1.5 and 1.8 arcmin 
respectively. 

Comparison between the X-ray images and radio (VLA- 
A observation of the H92a radio recombination line, Zhao et 
al. , 1997) and near-IR images (Wynn- Williams et al. , 1991), 
available in the literature, shows some clear similarities. As 
the radio and near-IR emission is clearly associated with the 
star-forming regions, the spatial coincidence implies that the 
soft X-ray emission has its origin in the starburst. We note 
that we do not see any evidence of X-ray emission from hot 
connecting the two galaxies, unlike the situation in some 
other interacting star-forming galaxies (e.g. the Antennae, 
Fabbiano et al. 1997). 



4 SPECTRAL ANALYSIS 

We fitted all the datasets simultaneously, thus covering an 
energy range from 0.1 to 8.0 keV. We rejected all ROSAT 
PSPC data points below 0.1 and above 2.0 keV; and used 
only the SIS and GIS data in the range 0.6 to 8.0 keV and 
from 0.7 to 7.0 keV respectively. In order to apply \ 2 statis- 
tics we re-binned the spectra so as to contain at least 20 
counts per bin. We used the software package XSPEC (vlO) 
to obtain the spectral fits. 

4.1 NGC3310 

We first fitted simple one-component models: a thermal 
bremsstrahlung model, a power-law model and finally a 
Raymond-Smith plasma model (Raymond & Smith 1977). 
These models were rejected at greater than the 99.9 per cent 
confidence level. The results of the spectral fits are given in 
table 5. 

We then fitted two component models: a two tem- 
perature Raymond-Smith plasma (ray-ray hereafter) and a 
Raymond-Smith plasma combined with a power-law (po- 
ray), following the results of Moran & Lehnert (1997) 
and Ptak et al. (1997) on the starburst galaxies M82 and 
NGC253. These models provided a much better fit than the 
simple models described above, at a confidence level over 90 
per cent for an addition of two free parameters (eg Bevington 
& Robinson 1992). We consider first the po-ray model. We 
obtained a photon index of T = 1.441°;^ and a temperature 
of kT = 0.8ll°;?2 keV with a reduced X 2 = 1-H5, where 
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Figure 2. X-ray contours from the HRI observation of NGC3310 overlaid on a POSS image. The contours are in levels of 2.2, 2.0, 1.7, 
1.33, 0.88, 0.44, 0.35, 0.26, 0.17 counts arcsec" 2 . 



Table 3. Sources in the HRI field of NGC3310 



RA (J2000) Dec (J2000) Count rate Flux Significance 

10 38 43.2 +53 31 07 2.0+0.27 x 10~ 3 2.3 x 10~ 14 15.5 

10 38 46.7 +53 30 38 2.8+0.26X 10~ 3 3.07 x 10~ 14 23.4 



Table 4. Sources in the HRI field of NGC3690 



RA (J2000) Dec (J2000) Count rate (cps) Flux (0.1-2. 5kcV) 



11 


28 


30.7 


+58 


33 


50 


10±1.2xl0~ 3 


1.18 x 10~ 13 


11 


28 


34.1 


+58 


33 


51 


15±1. 5xl0- 3 


1.72 x 10~ 13 


11 


28 


30.0 


+58 


34 


03 


3+0.6 xl0~ 3 


3.17 x 10~ 14 
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the abundance is fixed to the solar value and the absorb- 
ing column density to the Galactic value, Nh = 0.7 x 10 
cm~ , (Stark et al. 1992). Leaving the absorbing column 
density as a free parameter further improves the fit (reduced 
X 2 = 1.031); this is significant at a confidence level of > 99 
per cent for one additional parameter. In contrast, leaving 
the abundance as a free parameter does not improve the fit 
further at a statistically significant level. 

In the ray-ray case we do not obtain a significantly im- 
proved fit compared to the po-ray model (reduced \ 2 = 
1.022). Figure 4 shows the spectrum of NGC3310 with the 
double Raymond-Smith model. We obtained temperatures 
of 0.80lg;^ keV and lAMt 1 ^ keV for an absorbing hy- 
drogen column density of N H = 1.37t° f 2 x 10~ 20 cm" 2 . 
When we fixed the column density to the Galactic value 
(N H ~ 1.0 x 10 20 cm" 2 ) we obtained 0.84±°;^ and kT w 64 
keV. We also introduced two different absorbing column den- 
sities, for the soft and the hard components, following the 



analysis of M82 by Moran & Lehnert (1997) who showed 
that the hard X-ray component is obscured at the nucleus 
by Nh ~ 10 22 cm~ 2 . However, we do not find similar evi- 
dence here, since the hard component column has a value of 
N H ~ 1-77 x 10 20 cm~ 2 . Finally, we attempted to measure 
the element abundance in order to compare it with that for 
other star-forming galaxies (Sansom et al. 1996, Serlemit- 
sos et al. 1996), which appear to show a systematic trend 
towards sub-solar abundances. Unfortunately the data for 
our galaxies are not able to set any useful constraint on the 
abundances. 

We do not detect the FeK a line at 6.7 keV as would 
be expected for a hot thermal plasma. The low signal/noise 
data above ~ 5keV, does not allow us to set a useful upper 
limit on the line equivalent width. The lack of a strong FcK 
line could be explained by lower than solar abundances, as 
is probably the case in other well-studied starbursts such as 
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Double Thermal plasma 



channel energy (keV) 

Figure 4. The top panel shows ROSAT and ASCA spectra of 
NGC3310 with the best fit double Raymond Smith model. The 
bottom panel shows the ratio of the data points to the model. 

The SIS data are marked with dots while the PSPC data are 
marked with triangles. 





channel energy (keV) 

Figure 5. The top panel shows ROSAT and ASCA spectra of 
NGC3690 with the best fit double Raymond Smith model. The 
bottom panel shows the ratio of the data points to the model. 

The SIS data are marked with dots, while the PSPC data are 
marked with triangles. 



NGC253, (Ptak et al. 1997), and the Antennae, (Sansom et 
al. 1996). 

We note that a previous preliminary analysis of ASCA 
data by Serlemitsos et al. (1996), find a spewctral slope con- 
sistent with our value. However, they derive a considerably 
higher column density, but without the benefit of using the 
lower energy ROSAT data. 



4.2 NGC3690 

The spectral analysis follows the same procedures as for 
NGC3310. Simple one-component models such as a power- 
law and a Raymond-Smith plasma are similarly rejected 
(see table 6). Two-component models, such as a po-ray and 
a ray-ray model (with the hydrogen column density as a 



,+o.n 

-0.11 ' 

= 0.83+°' 0a 



kT = 0.83±°f 4 



, N H = 

+5.95 



0.03' ^2 — 10-312.44' 

20 cm" 2 respectively (the Galactic ab- 



free component) give T = 1.56 
2A±° 6 5 x 10 20 cm" 2 and kTi 

n h = i.6tg:2 

sorption is Nh = 1 x 10 20 cm -2 ). In the above model the 
abundance was fixed to solar. Also, as for NGC3310 when 
the abundance is a free parameter, it cannot be usefully con- 
strained. Figure 5 shows the ROSAT and ASCA spectrum 
of NGC3690 with the best fit double temperature Raymond- 
Smith thermal plasma model. We can clearly see from the 
spectrum that there is no evidence for an Fe line at ~ 6.7 
keV. The upper limit to the Fe line equivalent width is 
~ 860eV. We also tried to fit the models using a differ- 
ent absorbing column for the hard component. We obtain 
Nh = l.Ttg'g x 10 20 cm" 2 , which is comparable to the previ- 
ous best fit values, while the \ 2 is not significantly improved 
for one additional parameter. Thus there is no evidence for 
excess absorption in the hard component, in contrast to the 
results for M82 (Lehnert & Moran 1997). We also note that 
both the ray-ray and the po-ray models do not provide an 
adequate fit to the data as they can be rejected at over the 
99 per cent level of confidence using a \ 2 goodness of fit. 

We therefore attempted to fit more complicated mod- 
els, such as multi-temperature components, and a power- 
law distribution of temperatures. This is model 'cevmekl' in 
XSPEC, (see Done and Osborne 1997). This model gives a 



best fit of reduced \ 2 = 1-286, with a maximum temperature 
of ~ 57 keV, a power-law index of the temperature distri- 
bution of a = 0.221°;", and N H = 3.2ll°;^ x 10 20 cm" 2 . 
A combination of a power-law and the multi-temperature 
model gives a slightly better fit (reduced \ 2 — 1.211) with 



T = 1.521°;^, T max = 1.17l°; 2 g and a fixed to 1 - Finally, we 
considered using a non-equilibrium- ionization model (NEI) , 
such as that applied to hot gas resulting from supernova ex- 
plosions (Hughes and Singh 1994). This model assumes that 
the gas is instantaneously shock heated, and then applies 
corrections to the Raymond-Smith thermal plasma code to 
account for the non-equilibrium ionization fraction. An ad- 
ditional variable in the NEI model is nt, which is the product 
of the electron density and the time elapsed since the passage 
of the shock wave. In principle the results from fitting NEI 
models can be used to estimate the characteristic timescale 
over which the gas reaches ionization equilibrium. Clearly 
for a situation in which multiple supernova occur over an 
extended time period, it will not be possible to identify a 
unique timescale. However, NEI models have been applied 
to the optically thin plasma in galactic scale superwinds, 
such as observed in M82, Tsuru et al (1998). Unfortunately 
in the case of NGC3690, we do not have the benefit of mea- 
surements of the low energy lines of silicon and magnesium 
needed to quantify the importance of non-equilibrium ef- 
fects. For NGC3690 the NEI models give a best fit with 
a combination of a power-law and a NEI model. However 
the reduction in the x 2 is n °t statisticaly significant at a 
level above 90%, and so we conclude that our data does not 
justify further consideration of the NEI models. However, if 
new soft X-ray emission line data become available, it would 
be worthwhile returning to this question. 



5 DISCUSSION 

5.1 The soft X-ray emission 

The HRI images clearly show the extended X-ray emission 
in both NGC3310 and NGC3690/IC694. This supports a 
thermal origin for the soft emission, which may be due to 
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Table 5. Spectral fitting results for NGC 3310 



Parameter 



kT (KeV) 



Power-law Single Temperature Double Temperature Raymond Smith + 



N H (10 20 cm-' 20 ) 
X 2 I d.o.f. 
Fluxt (0.1-2.0keV) 
Fluxt (2.0-10.0kcV) 
Total Luminosity* (0.1-lOkcV) 



1 fiS+ 008 

o ic+0.63 
°- lo -0.52 

222.9/175 
0.95 
1.88 
1.3 



Raymond-Smith 

7 or + 1.24 



1 90+ ' 45 
239.3/175 

0.86 

2.12 
1.4 



Raymond Smith 

14 QS + 13 ' 52 



-I , 7 +0.50 
1 -°'-0.32 

168.7/165 
0.9 
2.10 
1.4 



Power-law 

,+0.09 
-0.12 



0.8T 



4+0.20 
-0.11 
,+0.68 
-0.40 



1.44] 
1.74 4 
170.2/165 

1.09 

2.10 
1.5 



■f All the fluxes arc in units of 10 ergs cm . |Thc luminosity is in units of 10 rmerg s 



a galactic scale super-wind eg Heckman et al. (1993). Such 
a phenomenon is observed in other star-forming galaxies, 
M82 (Strickland, Ponman & Stevens 1997), NGC1569 (Delia 
Ceca et al. 1996) and NGC4449 (Delia Ceca et al. 1997). 
According to this scenario, supernova in the starburst re- 
gion create a hot (~ 10 s K) bubble. The dense shell around 
the bubble fragments due to Rayleigh- Taylor instabilities 
and forms the optical emission line filaments and is prob- 
ably a source of soft X-ray emission. The optical data also 
support this model. In the case of NGC3690 the H a images 
show plumes and filaments extending out to several kpc (Ar- 
mus, Heckman & Miley 1990). Additional evidence comes 
from the extended synchrotron radio emission which may 
arise from the electrons produced by supernova remnants 
in the starburst regions, (Gehrz et al. 1983). In the case of 
NGC3310 the situation is more complicated as its low incli- 
nation angle (9 « 32°) makes detection of any outflow along 
the galaxy's minor axis very difficult. 

In order to further test the validity of the above model, 
we can estimate the expected X-ray luminosity of the super- 
wind and compare this with the observed luminosity of the 
soft Raymond- Smith component. First we relate the X-ray 
luminosity of the gas contained in a supperbubble with the 
basic properties of the starburst; its bolometric luminosity 
and age. Heckman et al. (1996) have calculated the total 
X-ray luminosity from a superbubble by integrating over its 
volume, and using the expressions for the gas density and the 
bubble radius derived by MacLow and McCray (1988). They 
find that L x ~ 5 x 10 39 L 33/35 n 17 / 35 tf /35 ergs" 1 , where 
Z/41 is the injection rate of mechanical energy into the bubble 
in units of 10 41 erg s _1 , n is the gas density and t-j is the age of 
the bubble in units of lOMyrs. Considering that the forma- 
tion of the bubble begins, as a result of strong winds from O 
stars following the starburst, we can assume that the age of 
the bubble is almost the same as the age of the starburst. Us- 
ing the relation dE/dt = 7 x 10 42 Lir i ii erg s _1 , (Heckman 
et al. 1996) which correlates the total infrared luminosity (in 
units of 10 11 L Q ) of a starburst with age greater than 10 7 yr 
with the kinetic energy deposition rate from supernovae and 
stellar winds, we can estimate the X-ray luminosity. Using 
the IRAS data from Soifer et al. (1987), the IR luminos- 
ity of NGC3310 is 5.5 x 10 10 L Q . Assuming a gas density of 
~ 1cm -3 and a typical starburst age of lOMyr we obtain 
an X-ray luminosity L x = 1.5 x 10 41 erg s _1 in the ROSAT 
band. This luminosity is larger than the measured luminos- 
ity in the same energy band of {Lx = 4.7x 10 40 erg s _1 ), but 
considering the large uncertainties in the model parameters 
involved, we conclude that the super-wind model remains a 



possible scenario for the origin of the soft X-ray emission. 
In NGC3690/IC694 the situation is more complex, as there 
are at least three distinct starburst components (Nakagawa 
et al. 1989) which are unresolved by IRAS. In order to ob- 
tain an estimate of the mechanical energy input for each 
component, we used the ground based mid-IR data (10-32 
microns) of Wynn- Williams et al. (1991). These observations 
have sufficient spatial resolution to separate the three com- 
ponents. We then make the reasonable assumption that the 
fractions of the total IRAS flux which arise from each of the 
three star-forming components, are similar to the fractions 
in the mid-IR band. Using the same formulae (Heckman et 
al. 1996) as in the analysis of NGC3310 above, and assum- 
ing a starburst age of 10 Myr and a density of n = 1 cm~ 3 , 
we predict L x = 1.3 x 10 42 , 1.1 x 10 42 and 3.2 x 10 41 erg s^ 1 
for the A,B and C components respectively (following nota- 
tion of Gehrz et al., 1983). Again the predicted luminosity 
is higher than the measured soft (0.1-2.5keV) X-ray lumi- 
nosity (4.7 x 10 41 crg s _1 ). This could possibly result from 
the fact that some fraction of the far-infrared emission mea- 
sured by IRAS, may arise from regions more extended than 
those producing the X-rays. Another possibility is that an 
additional source of heating may be present, perhaps due to 
stars of later spectral type than those producing the winds. 

The soft X-ray spectra of NGC3310 and NGC3690 are 
very similar. A Raymond-Smith component with a temper- 
ature of kT ~ 0.8 keV dominates the spectrum at soft ener- 
gies and could originate from the super-wind. This temper- 
ature is similar to that found for other star-forming galax- 
ies, eg M82 and NGC253 (Ptak et al. 1997), NGC1569 
(Delia Ceca et al. 1996), NGC4449 (Delia Ceca et al. 1997) 
NGC6240 and NGC2782 (Schulz et al. 1997). In the case 
of NGC3690 where the photon statistics are good, the SIS 
spectral fit yields a poor % 2 At- This could be suggestive of a 
hot tenuous gas component. Indeed, if the expanding gas has 
not had time to reach thermal equilibrium, its spectrum will 
not be well represented by the Raymond-Smith model. The 
use of a non-equilibrium code by Hughes et al. improved the 
fit, but again the % 2 is n °t accepted at the 98 per cent confi- 
dence level. This is not surprising as even a non-equilibrium 
model may be inadequate to fit an ensemble of supernova 
remnants occuring in several star-forming regions as is the 
case for NGC3690/IC694. 



5.2 The hard X-ray emission 

The spectrum at X-ray hard energies can be represented, 
in both galaxies, by a power-law (F ~ 1.4 — 1.6), or alter- 
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natively by a Raymond- Smith component with high tem- 
perature, (kT>10 keV). There is no strong evidence for the 
presence of a large amount of obscuration as is the case for 
NGC253 and M82 (Ptak et al. 1997). There are at least 
four possibilities for the origin of the hard X-ray emission: 
a) Inverse Compton scattering of the IR photons produced 
in the starburst by electrons from the numerous supernova 
remnants (eg Schaaf et al. 1989) b) emission from a low- 
luminosity AGN as is the case in NGC3628 (Yaqoob et al. 
1995) c) thermal emission from a very hot gas (T ~ 10 s K) 
and d) emission from X-ray binaries (Griffiths & Padovani 
1990). We discuss each of these possibility in turn. We con- 
sider first the possibility, that the hard X-ray emission arises 
from Inverse Compton (IC) scattering of the copious in- 
frared photons off of the relativistic electrons generated by 
the numerous supernovae. Support for this, especially in 
the case of NGC3310, comes from the similarity (within 
the errors) of the spectral index of the hard X-rays and 
the spectral index of the radio emission (a ra( j = 0.611q'q 3 
Niklas et al. 1997). Following Schaaf et al. (1989) we can 
estimate the X-ray luminosity from Inverse Compton scat- 
tering. Unfortunately, we can perform this calculation only 
for NGC3310, as there are no suitable radio data avail- 
able for NGC3690. From Vallee (1993) we have that the 
"minimum energy" magnetic field (see Longair 1992) is 
B = 0.47 x 10~ J Gauss. Then the minimum energy density 
is u = (7/3)(B 2 /8tt) erg cm" 3 = 2.05 x 10~ 12 erg cm" 3 . 
Following Schaaf et al. (1989) 

we have that Lie = (1/3)<ttRicLir^^- 72' 8 7i' 2 , where or, 
Ric, Lir and mc 2 , e e are the Thomson cross-section, the 
thickness of the disk, the far-IR luminosity, the rest mass of 
the electron and the energy density of the relativistic elec- 
trons, respectively. For a typical galactic disk Ric ~ lkpc 
(we cannot measure the actual thickness of the disk since 
NGC3310 is almost face on) and for a typical value of the 
low frequency cut-off v — 0.01GHz we have 71 = 150 for the 
lower limit of the Lorentz factor, and 72 = 10 for the maxi- 
mum Lorentz factor, in order to have IC emission at lOkeV 
(Schaaf et al. 1989). So for the Inverse Compton X-ray lumi- 
nosity we predict a value of 2.5 x 10 38 erg s" 1 , much lower 
than the detected hard X-ray luminosity of NGC3310, thus 
implying that IC can be only a minor component of the to- 
tal X-ray emission from this galaxy. We caution that this 
result is quiet uncertain as the calculation depends on pa- 
rameters like the volume of the source and the thickness of 
the disk which are poorly known. Another possible origin of 
the hard X-ray emission is the presence of a low luminosity 
AGN, although there is no evidence for non-starforming nu- 
clear activity from diagnostic emission line ratio diagrams 
based on optical and near infrared spectra. 

However, the hard X-ray power-law spectrum is never- 
theless consistent with the presence of an active nucleus. In 
order to test this possibility further we use the Lx/Lu a 
relation from Elvis et al. (1984), where Ln a is the lu- 
minosity of the broad component of the H a line and L x 
is the hard X-ray luminosity (2-10 keV). Lx/Ln a — 40, 
for low- luminosity AGN and thus we estimate f(H a ) ~ 
5.3 x 10~ 14 erg s" 1 cm" 2 and 2.7 x 10~ 14 erg s" 1 cm" 2 for 
NGC3310 and NGC3690 respectively. A broad component 
of this strength is easily detectable, but is not seen (Ho 
ct al. 1997). The absorbing columns found from the X-ray 
spectral fitting (ray-po models in tables 5 and 6) imply low 



extinctions, Av ~ 0.1 mag for both galaxies (Bohlin et al. 
1978). Hence a BLR reddended by this amount would still 
be observed. 

The next possibly is the presence of a very hot thermal 
component. A Raymond-Smith model with kT>10 keV pro- 
vides a good fit to the data. However, the strong FeK line 
at 6.7 keV which should accompany the thermal emission is 
not observed. This could be attributed to a low metallicity, 
but unfortunately we cannot check this as the abundances 
are not well constrained using the present X-ray data. Al- 
ternatively, the lack of a FeK line could be explained by a 
low contribution of type la supernovae to the enrichment of 
the Interstellar Medium (Sansom et al. 1996). 

Finally we consider high mass X-ray binaries as a possi- 
ble origin of the hard X-ray component. High mass binaries 
will form as a consequence of the starburst, and indeed many 
such systems have been identified in nearby star-forming 
galaxies (Read et al. 1997, Fabbiano 1995 and references 
therein). Assuming that a typical X-ray luminosity of these 
systems is io 37-38 erg s _1 , we can estimate the total num- 
ber of binaries. As the X-ray luminosity of the hard compo- 
nent is L x ~ 5 x 10 41 and L x ~ 10 41 erg s" 1 for NGC3690 
and NGC3310 respectively, we estimate a range of between 
5,000-50,000 and 1,000-10,000 X-ray binaries for the two 
galaxies. Now we can compare this with the number of ion- 
izing OB stars determined from the integrated far-infrared 
luminosity. Making the assumption that it is mostly these 
stars that heat the dust, which then reradiates in the mid- 
far infrared, we estimate ~ 2x 10 5 and ~ 3.5 x 10 6 OB 
stars for NGC3310 and NGC3690, respectively. If, follow- 
ing Fabbiano et al. (1992), 0.2 percent of these are massive 
X-ray binaries, then there are 400, and 7000 such systems, 
respectively. For the upper range of binary luminosities the 
predicted X-ray luminosity is comparable to that observed 
in both cases. Indeed, if the hard X-ray emission arises from 
binaries with low-metallicity and thus high X-ray luminosity 
(L x ~ 20 erg s _1 ) like those observed in the Magellanic 
Clouds (van Paradijs & McClintock 1995), then they could 
easily produce the observed luminosity. One potential prob- 
lem is that some point sources observed in nearby galaxies 
with the ROSAT PSPC by Read et al. (1997), appear to 
have soft spectra (kT ~ 2 keV), whilst the the high tem- 
peratures inferred for NGC3310 and NGC3690 are closer to 
those predicted for high mass X-ray binaries (see Nagase 
1989). 



6 CONCLUSIONS 

We have modelled and interpreted the combined ASCA and 
ROSAT X-ray spectra from 0,1-10 keV, for the star-forming 
galaxies NGC3310 and NGC3690. These two galaxies are 
amongst the most X-ray luminous [L x ~ 10 41 ~ 42 ) erg s" 1 
starbursts in the local Universe. In addition, ROSAT HRI 
images show that the emission from NGC3310 is extended 
out to at least ~ 3 kpc. The soft X-ray emission from 
NGC3690 comes from at least three spatially resolved re- 
gions. The limited spatial resolution of the ASCA data does 
not allow us to place useful limits on the extent of their hard 
X-ray components. The X-ray spectrum of NGC3310 can be 
described by two components: at soft energies a Raymond- 
Smith component (kT~0.8 keV) which probably originates 
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from a super-wind. The predicted soft X-ray emission, on 
the basis of the supernovae mechanical energy deposited to 
the interstellar medium, is comparable to that observed. At 
harder energies we can fit either a Raymond-Smith compo- 
nent (kT~17 keV) or a power-law T ~ 1.4. 

The results for NGC3690 are similar. There are at least 
two components in its spectrum: a soft Raymond-Smith 
plasma (kT~0.8 keV) and a harder component which can 
be represented equally well by a Raymond-Smith compo- 
nent (kT~10 keV) or a hard power-law V ~ 1.6. However, 
the best-fit model for NGC3690 is rejected at the 98 per cent 
confidence level implying that more complicated models are 
necessary. Hence we considered a non-equilibrium ionization 
model (Hughes et al. 1994) and a multi-temperature ther- 
mal model (Done & Osborne 1996). The above models do 
improve the fit, but only at the ~ 90 per cent confidence 
level. Although non-equilibrium models may be important, 
additional information such as X-ray emission line ratios are 
required in order to make further quantitative progress. 

In neither NGC3310 and NGC3690 do we find evidence 
for significant absorption columns above the Galactic values, 
in either the soft or the hard X-ray component. 

The nature of the hard X-ray emission is still uncertain, 
and we considered various possibilities. Of these we conclude 
that Inverse Compton scattering by high energy radio elec- 
trons and infrared photons, can probably provide only a mi- 
nor contribution. Although the presence of an AGN is a pos- 
sibility, this is not supported by data at other wavelengths, 
and the modest columns argue against an obscured nucleus. 
A very hot thermal component may be present, (~ 10 s K). 
We do not detect the FeK line at 6.7 keV, possibly because 
the the abundances are sub-solar. To test this possibility 
further requires information on element abundances, and an 
extension of the observed spectra to higher energies. Finally, 
estimates of the number of X-ray binaries, based on the to- 
tal mid-IR emission, suggests that they may well be able to 
account for the hard X-ray emission. 

The X-ray results on these two galaxies are similar 
to those obtained for other dwarf star-forming galaxies 
(eg NGC1569, NGC4449) and the archetypal star-forming 
galaxies M82 and NGC253, although there is a range in rel- 
ative contributions from the various X-ray components. This 
implys that the same general emission mechanisms apply in 
star-forming galaxies over three decades of luminosity. 
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